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ABSTRACT 

A  preliminary  simulation  model  of  the  search  mode  of  a  phased 
array  radar  has  been  developed.  The  basic  components  are  a  search 
strategy  for  the  radar,  an  evasion  policy  for  the  target,  and  a 
simple  radar  environment.  At  present  there  is  no  attempt  at  adap¬ 
tive  modification  of  either  component.  A  particular  search  strategy 
has  been  chosen.  There  is  provision  in  the  simulator  to  modify  the 

f 

procedure,  however.  The  evasion  policy  of  the  target  is  one  of 
changing  course  at  random  times  with  random  deviations. 

It  is  planned  to  develop  the  tracking  package  next.  The 
approach  will  be  through  a  Kalman-Bucy  filter. 


PROBLEM  STATUS 

This  is  an  interim  report;  work  on  this  problem  is  continuing. 


AUTHORIZATION 

NRL  Problem  No.  B01-07 
NAVSHIPSYSCOM  Project  SF  001-02-05-6284 
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RESOURCE  MANAGEMENT  OF  FLEXIBLE  SENSORS 
I.  PHASED  ARRAY  RADAR  SIMULATOR 


INTRODUCTION 

A  computer  program  has  been  developed  for  the  CDC  3800  computer 
to  simulate  a  radar  detection  system  and  target.  The  radar  detection 
system  is  located,  by  assumption,  on  a  ship  at  a  specified  position 
and  orientation  on  the  earth  and  will  be  referred  to  as  the  obser¬ 
ver.  This  system  scans  particular  ranges,  elevations,  and  azimuths 
according  to  a  predetermined  strategy.  The  pulse  rate  and  power 
transmitted  by  the  system  at  any  given  time  are  also  specified  by 
the  strategy.  The  target  is  assumed  airborne  and  flying  at  con¬ 
stant  altitude  and  speed,  heading  by  way  of  an  evasion  course 
towards  the  observer.  The  observer  and  target  are  considered 
adversaries,  and  the  present  goal  is  for  the  observer  to  detect  the 
target  before  the  target  is  close  enough  to  destroy  the  observer. 

DESCRIPTION  OF  SIMULATOR 

Figure  1  is  a  general  block  diagram  of  the  simulator.  It  will 
be  described  briefly  in  the  remainder  of  this  section.  The  portions 
of  the  block  diagram  indicated  by  asterisks  are  each  detailed  pro¬ 
cedures  and  will  be  described  in  subsequent  sections. 

For  simplicity,  the  initial  model  of  the  radar  beam  is  a  conical 
surface  with  apex  at  the  observer,  bounded  above  and  below  by  an 
upper  and  lower  elevation  (fi  and  6„)>  and  bounded  to  t:he  sides  by 

U  Xj 

an  upper  and  lower  azimuth  (vy  and  v^),  where  6^  ■  vy-  ■  1.7°. 
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This  model  will  be  increased  in  complexity  after  other  more  primary 
goals  are  attained,  in  particular  the  incorporation  of  target  track¬ 
ing  into  the  simulator. 

The  target  is  considered  a  point  in  space,  with  range,  elevation, 
and  azimuth  coordinates  at  time  t  expressed  as  RT(t),  6T(t),  and 
v  (t)  respectively.  The  sequence  of  operations  is  roughly  as  follows, 
The  main  program  calls  on  a  strategy  subroutine,  and  the  latter 
yields  the  output  parameters  of  time  t,  beam  position  an<* 

v  )  at  time  t,  and  the  relative  power  radiated  (mode  of  operation). 

The  main  program  then  summons  the  target  simulator  subroutine  to 
obtain  the  target  coordinates  RT(t),  6T(t),  vT(t)  at  the  time  t 
specified  by  the  strategy.  A  test  is  then  performed  to  see  if 
6t  >  65°,  if  this  is  the  case,  then  the  target  is  considered 
within  destruct  distance  of  the  observer  and,  accordingly,  the 
observer  is  assumed  destroyed  (it  should  be  remarked  that  flT  >  65° 
is  an  arbitrary  and  perhaps  somewhat  inaccurate  criterion  for 
destruction,  but  it  is  being  used  temporarily  for  expediency). 

This  terminates  the  program.  If,  on  the  other  hand,  £  65°,  the 
program  continues  and  next  determines  whether  the  target  is  inside 
or  outside  of  the  radar  beam.  It  is  inside  if  both  inequalities 
8^  s  5?  s  flu  an<*  s  vt  S  vu  are  8at*ifie{*>  ^  this  *8  case, 
then  a  probability  of  detection  is  computed  based  on  the  range  of 
the  target  and  its  position  in  the  beam,  and  on  the  mode  of  the  beam, 
A  random  number  between  0  and  1  is  chosen  and  compared  with  the 
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probability  of  detection:  if  the  former  is  no  greater  in  value  the 
target  is  assumed  discovered,  and  the  program  terminates;  otherwise, 
control  cycles  back  to  the  calling  of  the  strategy  subroutine  to 
obtain  the  next  set  of  parameters.  Control  also  cycles  back  to  this 
point  if  the  target  la  outside  the  beam. 

The  following  sections  will  describe  in  more  detail  the  strategy 
and  target  simulator  subroutines,  and  the  computation  of  the  proba¬ 
bility  of  detection.  The  results  of  a  series  of  computer  tests  of 
the  general  simulator  can  be  seen  in  the  Appendix. 

TARGET  SIMULATOR 

The  method  of  operation  of  the  target  simulator  subroutine  will 
now  be  described.  This  Involves  a  set  of  assumptions  about  the 
earth  and  target,  a  method  of  choosing  the  type  of  motion  of  the 
target,  including  evasive  action,  and  finally  the  derivation  of  the 
equations  of  target  position  at  some  particular  time  t. 

Assumptions  about  Target  and  Environment 

The  earth  is  assumed  to  be  a  sphere  of  radius  R_  ■  3963.34  stat¬ 
ute  miles  (the  term  "mile(s)"  will  subsequently  mean  "statute  mile(s)" 
unless  otherwise  specified),  and  the  observer  is  postulated  to  be  at 
the  approximate  location  of  Washington,  D.  C.:  colatitude  52°,  east 
longitude  283°, 

Because  the  target  remains  at  a  constant  height  h  above  the 
earth's  surface,  its  motion  is  restricted  to  the  surface  of  a  sphere 
having  radius  R^  +  h  and  concentric  with  the  earth.  If  we  radially 
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project  the  observer's  coordinates  onto  this  sphere,  hereafter  referred 
to  as  the  target  sphere ,  we  can  analyze  the  relation  between  the  target 
and  projected  observer  usins  spherical  t-r-f  trnnnrnAi-ry  n  stipulated 
that  the  target  first  appears  on  the  north  horizon  of  the  observer;  it 
thus  has  a  colatitude  of  less  than  52°  and  an  east  longitude  of  283°. 
The  evasion  course  of  the  target  consists  of  a  series  of  connected 
great-circle  paths  or  segments;  two  consecutive  segments  will  be  said 
to  intersect  at  a  pivot  point .  While  the  direction  of  the  target 
changes  abruptly  at  these  pivot  points  its  speed  is  constant.  At  any 
particular  pivot  point  the  angular  deviation  of  the  next  great-circle 
path  from  a  direct  great-circle  path  to  the  projected  observer  is 
chosen  at  random  according  to  a  Gaussian  distribution.  The  time  that 
the  target  adheres  to  any  particular  great-circle  path  (i.e.,  the 
time  between  pivot  points)  is  also  chosen  at  random  according  to  the 
waiting  time  distribution  of  a  Poisson  process. 

Equations  of  Target  Position 

The  following  symbols  are  used  in  developing  the  equations  of 
position  of  the  target  at  some  time  t: 

tQ,  t^,  t2>  t^,  ...  are  the  times,  in  sequence,  at  which 

the  target  arrives  at  each  pivot  point.  At  time  t  the  target  is  at 
the  north  horizon  of  the  observer;  t^  is  the  time  of  arrival  at  the 
first  pivot  point  thereafter;  t^  the  time  of  arrival  at  the  next  pivot 
point,  and  so  forth.  Clearly  t0  <  <  C2  <  t3  <  •••  <  t^  <  ...  ■ 
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i  i  th 

The  point  (0q,0o)  gives  the  coordinates  of  the  target  at  the  i 

pivot  point  (occurring  at  time  t^,  where  0*  and  0*  are  the  target 
colatitude  and  east  longitude,  respectively.  Thus  the  target  path 
is  a  series  of  great-circle  arcs  joining  the  points  (0*  ^,0^  S  and 


(eJV),  i  -  i,  2,  ...  . 


0^  is  the  deviation  the  target  path  makes  at  the  point  (0*,0*) 
from  a  direct  great-circle  path  Joining  (0*,0^)  and  the  projected 
observer.  If  the  deviation  is  counterclockwise  (as  viewed  from 
outside  the  target  sphere)  a  0,  if  clockwise  then  <  0. 

The  point  (0  ,0  )  gives  the  coordinates,  colatitude  and  east 

.D  D 

longitude  respectively,  of  the  observer  (and  thus  also  the  pro¬ 
jected  observer).  Hence  0  ■  52°  and  0_  ■  283°. 

B  B 

h  is  the  height  of  the  target  in  miles. 

VT  is  the  speed  of  the  target  in  mph. 

a)  is  the  angular  velocity  of  the  target  in  rad/sec.  We  have 
(  3600 


that 


VT 

vh 


The  point  (0,0)  ■  (0(t),0(t))  gives  the  colatitude  and  east  lon¬ 
gitude  of  the  target  at  time  t. 

The  range,  elevation,  and  azimuth  of  the  target  at  time  t  rela¬ 
tive  to  the  observer  are  denoted  respectively  by  RT(t),  fiT(t),  and 
vT(t).  The  argument  "(t)"  will  sometimes  be  omitted  for  convenience. 

Suppose  at  time  t  the  target  has  travelled  a  slight  distance 
beyond  its  most  recent  pivot  point  which  it  reached  at 
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time  tj.  Thus  t^  <  t  &  t^+^,  and  c^e  great-circle  path  the  target  is 
now  traversing  makes  an  angle  a ^  with  the  great-circle  arc  joining 
(0^)  and  (0B,0B).  Figure  2  is  the  representation  of  this  situa¬ 
tion  on  the  surface  of  the  target  sphere.  We  can  calculate  the 
target  coordinates  0(t),  0(t)  using  the  known  values  0  ,  0  ,  0*,  0*, 
uo ,  t,  t,,  and  The  following  procedure  accomplishes  this: 

Step  1:  Calculate  p  from  the  relation 

p  ■  cos-1[cos  0B  cos  0*  +  sin  0_  sin  01  cos  (0*  -  0n)]. 

D  O  D  O  O 

Step  2a:  Calculate  cos  0  from  the  relation 

„„„  a  C0S  ‘  cos  ^  cos 
COS  p  m  _ B _ O 

sin  p  sin  0* 

If  |  cos  0 1  £  0.999,  0  is  computed  as  the  arccosine  of  cos  0,  noting 
that  0  s  p  s  180°,  and  we  proceed  to  step  3;  otherwise  we  perform 

Step  2b:  If  |  cos  0|  >  0.999,  calculate  0  from  the  relation 


(1)  0  -  90°  +  [sgn (cos  0)]  atnlf*in  eBl8ln®o  "  0B^V  90°  , 

\  sin  p  ' 

in  order  to  reduce  numerical  error.  The  principal  arcsine  is  used,  and 
the  "sgn"  function  is  defined  as  follows: 

sgn  X  s  +1  if  X  2  0, 
s  -1  if  X  <  0, 

where  X  is  any  real  number.  If  relation  (1)  is  used,  0  satisfies 
the  inequality  OS  p  s  180°. 

Step  3:  Compute  y  from  the  formula 

y  -  0  +  *  { sgn  [sin  (0^  -  0fi)3}  .  0sy<  360°, 
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Step  4:  Compute  0(t)  from  the  formula 

”11  i 

9(t)  ■  cos  [cos  0^  cos  ^(t-t^)  +  sin  0q  sin  uKt-t^cos  y] 

St  >  5a:  Compute  cos  l|)  from  the  relation 

cos  u>(t-t.)  -  cos  0(t)  cos  0* 

cos  Ip  ■  - ; - - - 2 

sin  0  (t)  sin  01 
o 

If  Icos  Ip  i  £  0.999,  tp  is  computed  as  the  arccosine  of  cos  Ip,  noting 
that  0  £  ip  £  180°,  and  we  proceed  to  step  6;  otherwise  we  perform 
Step  5b:  If  Icos  Ip  l>  0.999,  calculate  Ip  from  the  relation 


Ip  ■  90°  +  [sgn(cos  $)] 


■  in'1/"8lnw(t’ti>|BlnylN'-  90 


o 


sin  0(t) 


■) 


where,  as  in  eq.  (1),  the  arcsine  is  taken  in  the  principal  region. 

Again  we  have  0  £  Ip  £  180°,  We  assume  also  that  the  pivot  points 

occur  sufficiently  often  to  insure  that  lu(t-t^)  <  rr radians. 

Step  6:  Compute  0(t)  from  the  relation 

0(t)  -  0*  -  ip  •  [sgn(siny)]  •  {sgn[sin(0*  -  0fi) ]}  , 

0(t)  should  satisfy  0  £  0(t)  <  360°. 

i+1 

We  note  that  if  t  »  t^+^  we  have  computed  0Q  ■  and 

i+1 

0Q  “0(ti+1),  ^e  coordinates  of  the  (i+l)at  pivot  point;  these 
coordinates  are  stored  for  possible  use  in  later  calcualtiona. 

We  digress  at  this  point  to  discuss  how  the  times  between  pivot 
points,  At  ■  t^  -  t^  (i«0, 1,2, . . . ),  and  the  deviations  0^ 
(i"0,l,2, .  . .) ,  are  determined.  To  choose  a  value  of  A^  for  any  i, 
we  must  first  have  access  tc  a  preselected  value  of  the  mean  free 
time  of  traverse  between  pivot  points,  denoted  by  tf,  which  remains 
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constant  for  the  entire  simulation.  A  random  number  r  satisfying 
0  <  r  <  1  is  then  chosen  from  a  uniform  distribution.  The  value  of 
At^  is  then  computed  as  At^  *  t^in  (  ).  The  value  of  tr  cur¬ 

rently  being  used  is  20  sec.,  which  guarantees  with  95%  probability 
that  At^  £  60  sec.  The  deviations  are  chosen  from  a  Gaussian 
distribution  having  a  preselected  standard  deviation  a ,  which 
remains  constant  for  the  entire  simulation.  To  date  we  have  per¬ 
formed  simulations  using  a a  ■  15°  and  a a  ■  30°,  using  the  CDC  COOP 
subroutine  G5  WISC  RANSS  to  generate  the  a^. 

We  now  calculate  the  values  R^t),  6 ^ ( t ) ,  and  vT(t)  at  time  t 
using  the  parameters  we  have  previously  specified  or  calculated. 
Figure  3  shows  a  meridian  plane  of  the  earth,  the  particular  plane 
being  the  one  containing  the  target.  The  following  procedure  then 
follows  Immediately: 

Step  1:  Calculate  £  from  the  relation  (cf .  Fig,  2) 

£  ■  cos  1  [cos  cos  0(t)  +  ain  0R  sin  0(t)  cos  (0(t)-0R)]. 

Step  2:  Calculate  R^  from  the  relation 


Rt  -  Vh2  +  2RE(RE+h)(l  -  cos  £). 


Step  3:  Observing  6T  ■  )?  -  90  ,  we  obtain  an  equation  for  com¬ 


puting  fiT, 


fiT  *  «in 


h(2RE+h)  -  r£ 


90°  s6ts  90°, 


8 


CONFIDENTIAL 


CONFIDENTIAL 


although  cases  where  5^,  <  0°  mean  that  the  target  is  beneath  the 
horizon. 

C  A.  f/t  T?  4  /V  H  f.A 

^  ~r  *•  *kw',“  w*-  **  --“C3  ~  —  *--£3'  -  I  ••  —  **«*-«•*-** 


coa  .  cos  8  (t)  -  cos  9b  cos  £  ' 
sin  0R  sin  £ 

If  Icos  vT!  £  0.999,  we  define  v'  to  be  the  principal  arccosine  of 

cos  v  so  that 
T 

■  cos  ^  (cos  vT),  0  s  v;  s  180°. 


We  then  compute  from  the  relation 

vT  -  180°  +  {sgn[sin(0<t)  -  0^)]}  •  (v '  -  180°) . 

If  Icos  v  |  >  0.999,  then  vT  is  computed  using  the  formula 


vT-  90  +  [sgn (cos  vT)] 


|ln«l /sin  0(t)  sin  (0(t)  -  0^N 
V  sin  £  J 


90°  » 


where  the  arcsine  is  taken  to  be  in  the  principal  region  [-90°, 90°]. 
vT  ia  normalized,  if  necessary,  so  that  0  £  vT  <  360°. 

STRATEGY 


We  subsequently  describe  the  present  strategy,  keeping  in  mind 
that  all  aspects  of  it  are  flexible  and  can  be  completely  modified 
for  future  experimentation.  The  present  strategy  was  based  on  con¬ 
ventional  procedure  and  was  employed  to  hasten  the  initial  testing 
of  the  simulator. 
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Antenna  Structure  and  Search  Modes 

The  present  model  of  the  observer  la  a  ship  with  four  antennas, 
each  with  a  90°soan  in asimuth.  mounted  In  tun  unite  «ith  »n  «ntenn« 
pair  in  each  unit  in  such  a  way  that  the  units  span  non-overlapping 
180°  sectors  in  asimuth.  Figure  4  indicates  the  positions  of  these 
antennas  on  the  ship.  The  antenna  pair  shown  on  the  right-hand 
part  of  the  diagram  spans  an  asimuth  range  from  0°  to  180°.  This 
pair  shall  be  subsequently  considered  as  one  antenna  for  simplicity; 
in  fact,  it  will  be  denoted  aa  antenna  #1.  Similarly,  the  remaining 
pair  spans  an  asimuth  range  from  180°  to  360°  and  will  be  denoted 
aa  antenna  #2. 

Three  search  modes  are  used  in  the  strategy,  and  these  shall 
be  designated  as  modes  I,  II,  and  III.  The  region  of  space  that 
is  searched  by  each  mode  will  also  be  designated  by  the  Roman 
numeral  of  the  mode.  Thus  we  can  refer  to  regions  and  modes  inter¬ 
changeably.  We  stipulate  that  the  power  transmitted  for  searching 
in  modes  I  and  III  is  the  same,  and  equal  to  one-fourth  that  trana- 
mitted  in  mode  II.  This  condition,  of  course,  will  penalise  the 
radar's  performance.  It  may  arise,  however,  under  certain  types 
of  energy  management,  e.g.,  purposely  transmitting  several  main 
antenna  lobes  simultaneously.  Future  versions  of  the  simulator 
will  provide  for  varying  degrees  of  such  management. 

We  define  R*,  R^1,  and  rJ11,  corresponding  to  the  respective  modes, 
as  those  ranges  at  which  the  target,  when  in  the  center  of  the  beam, 
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can  be  detected  with  a  507o  (0.5)  probability.  We  specify  that  R^,  *  150 
nautical  miles  (n.m.).  The  above  conditions  determine  R*  and  R^11 


\  Icjuc  ly 


4  r.Vi  4  r 

»V  Vi.  ViiV  i.  MVU4.  AV14|  »“■*■-**  —  ~ 


Its  fundamental  form  is 


ptGtV 

(4W*‘)2 


where  P  ■  echo  power  received  at  the  radar, 

P  ■  power  transmitted  by  the  radar, 

Gfc  ■  transmitting  gain  of  the  antenna. 

A  ■  effective  capture  area  of  the  receiving  antenna, 

a  ■  radar  cross  section  of  the  target, 

and  R^,  -  range  of  the  target. 

Let  superscripts  I,  II,  and  III  refer  to  quantities  relating  to 

modes  (or  regions)  I,  II,  and  III,  respectively.  The  variables  Af 

and  a  are  independent  of  the  mode  number,  and  the  antenna  gain  Gt 

is  kept  independent  by  considering  the  target  to  be  in  the  center 

1  pj 

of  the  beam  in  all  cases.  Thus  P'  ■  K  t  ,  where  j  ■  I,  II,  or 

(4>4 

III,  and  K  is  a  lumped  constant  independent  of  J.  Since  equal 
detection  probabilities  mean  that  equal  echo  power  is  received 
at  the  radar,  we  have  P*  -  P*1  -  P*11.  We  then  get  for  the  radar 


equation 


pI  pH 


,  or  -  R^1 
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Since  P*1  ■  4P*  and  R*1  ■  150  n.m.,  it  follows  from  (2)  that 

R*  ■  106.1  n.m,  ■  R*^,  the  last  equality  occurring  because 

pi  .  piii. 

Figure  5  depicts  these  regions  In  polar  coordinates,  with  range 
as  the  modulus  and  elevation  the  argument.  The  minimum  and  maximum 
elevations  searched  in  each  region  are  Indicated,  as  are  the  ranges 
■}(-  R*11)  snd  R^1.  The  diagram  also  shows  the  frequency  with  which 
each  region  is  searched.  For  any  particular  elevation  within  the 
bounds  searched  in  any  of  the  modes,  all  aslmuths  from  0°  to  360° 
are  scanned.  Region  I  is  completely  searched  every  second,  regions 
II  and  III  every  12  and  6  seconds  respectively.  Elevations  from 
0°  to  1,7°  (one  beamwldth)  are  scanned  in  region  I,  0°  to  4.5° 
in  region  II,  and  1.7°  to  65°  in  region  III. 

The  pulse  rate  used  in  each  region  is  a  function  of  the  maxi¬ 
mum  range  to  be  searched  in  that  region.  This  dependency  is 
brought  about  by  the  conflicting  goals  of  having  as  high  a  pulse 
rate  as  possible  to  minimize  search  time,  yet  slow  enough  so  that 
the  return  pulse  will  be  received  before  the  next  transmitted  pulse 
occurs  and  thus  not  create  ambiguity  as  to  the  source  of  the  return 
pulse  and  the  value  of  target  range.  The  maximum  unambiguous  pulse 
rate  f  is  related  to  the  maximum 


range  by  the  formula 


2rmax 
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where  C  is  the  velocity  of  light.  We  use  as  the  value  of  R^^  the 
ranges  R*,  R*1,  and  R*11,  at  which  the  probability  of  detection  in 

T  T  TT 

each  respective  mode  is  0.5.  It  follows  that  f‘  *  f“*  *  761  pulses/ 

R  R 

second,  and  f^  ■  538  pulaea/seeond.  Accordingly  s  pulse  rate  of 
537 /sec.  was  used  for  region  II  and,  through  an  oversight,  a  pulse 
rate  of  1253/sec.  was  used  for  regions  I  and  III  (this  is  not 
expected  to  influence  significantly  the  results  of  our  simulator 
teats,  because  we  are  not  determining  target  range  at  this  stage  of 
the  project;  however,  this  error  will  be  corrected  prior  to  subse¬ 
quent  tests).  These  rates  suggest  we  create  a  new  unit  of  time, 
denoted  as  a  ana if ,  such  that  1  second  *  3,759  snarfs.  Then  for 
region  II,  with  a  maximum  pulse  rate  of  537/sec.,  we  must  have  at 
least  7  snarfs  between  pulses.  For  regions  I  and  III,  with  a  maxi¬ 
mum  pulse  rate  of  1253/sec.,  there  must  be  at  least  3  snarfs  between 
pulses.  Our  strategy  operates  at  these  maximum  pulse  rates. 
Representation  of  Beam  Positions 

The  beam  positions  used  in  the  present  strategy  are  considered 
discrete;  that  is,  the  beam  does  not  move  continuously  through  a 
region,  but  Instead  scans  a  region  by  Incremental  steps  in  position, 
either  in  elevation,  asimuth,  or  both,  Thus  at  most  a  finite  number 
of  possible  beam  positions  can  be  assumed  by  the  radar.  We  mention, 
however,  that  the  particular  choice  of  discrete  beam  postlons  used 
here  is  temporary  and  probably  will  be  modified  in  the  future  to  pro¬ 
vide  a  more  tightly  packed  scan  of  each  region. 
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The  azimuth  incremental  positions  used  are  independent  of  the  region. 
For  any  particular  azimuth  scanned  by  antenna  #1,  an  azimuth  180  degrees 
greater  is  being  scanned  by  antenna  #2.  Thun  if  a t  any  instant  of  time 
the  beam  center  of  antenna  #1  is  pointing  at  an  azimuth  v,  the  beam 
center  of  antenna  #2  is  pointing  at  azimuth  v  +  180°.  Each  azimuth 
position  can  be  represented  by  the  format  [v^ , vy] ,  where  and 
are  expressed  in  degrees.  We  use  this  in  spite  of  the  redundancy  aris¬ 
ing  from  the  relation  vu  “  +  1-7°.  There  are  106  positions  in 

azimuth  that  each  antenna  can  assume  under  the  present  strategy,  cover¬ 
ing,  in  slightly  overlapping  segments,  the  span  of  0°  to  180°.  These 
are  the  positions 

[1.6981(n-l) ,  0.0019  +  1.6981n],  n  -  1,  ....  106. 

This  is  the  minimum  number  of  positions  necessary  for  an  antenna  with 
1.7°  beaimrldth  to  span  180°  in  azimuth.  Corresponding  to  the  positions 
of  antenna  #1,  antenna  # 2  assumes  the  positions 

[180  +  l,6981(n-l),  180.0019  +  1.6981n],  n  ■  1,  ....  106. 

Thus  an  azimuth  position  is  described  by  the  antenna  number  and  the 
integer  n. 

The  elevation  incremental  positions  used  are  dependent  on  the 
region.  Similar  to  the  azimuth  format,  we  represent  the  elevation 
by  [6., 6  ],  where  ft.  and  5  are  expressed  in  degrees  and  8„  *  6.+  1.7°. 

In  region  I  there  is  only  one  possible  elevation  position,  and  that 
is  [0,1.7],  In  region  II  there  are  three  possible  elevation  positions: 

[1.4(m-l) ,  0.3  +  1.4m],  for  m  -  1,  2,  3. 
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In  region  III  there  are  38  possible  elevation  positions,  and  they  are 
[1.6658(m-l) ,  0.0342  +  1.6658m],  form-  2,  3,  ....  39. 

In  each  of  regions  li  and  til  tne  positions  are  slightly  overlapping 
and  scan  the  total  interval  using  the  minimum  number  of  positions 
necessary  for  a  beamwidth  of  1.7°.  Antennas #1  and  #2,  although 
180°  apart  in  asimuth,  are  always  searching  at  the  same  elevation. 

The  preceding  discussion  indicates  that  the  beam  position  and 
power  transmitted  can  be  completely  specified  by  an  ordered  quad¬ 
ruple  of  integers;  (1)  the  region  no.,  (2)  the  antenna  no.,  (3)  m, 
and  (4)  n  (for  region  I,  m  is  set  to  1). 

Sequence  of  Beam  positions  and  Modes  for  Current  Strategy 

The  search  sequence  used  in  the  present  strategy  has  a  period  of 
12  seconds.  At  the  beginning  of  the  0th  second,  the  scan  of  region 

I  begins  at  pulse  rate  f*  and  continues  to  completion.  7  snarfs 

R 

later,  the  scan  of  region  II  begins  at  pulse  rate  f**  and  continues 

until  one-twelth  of  this  region  has  been  searched.  3  snarfs  later, 

ITT  I 

the  scan  of  region  III  begins  at  pulse  rate  fR  ■  fR  and  continues 
until  one-sixth  of  this  region  has  been  searched.  At  this  point 
about  0.7  sec.  have  elapsed,  and  no  further  searching  1b  done  until 
the  beginning  of  the  1st  second.  Then  the  scan  of  region  I  begins 
again  and  is  carried  to  completion  as  before.  7  snarfs  later  the 
scan  of  region  II  resumes  at  the  point  at  which  it  was  interrupted 
during  the  0th  second,  and  another  twelth  of  this  region  is  searched. 
3  snarfs  thereafter,  the  scan  of  region  III  resumes  at  the  point  of 
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interruption  during  the  O'-*1  second,  and  another  sixth  of  this  region 
scsrchcd.  Tiisn  a«-aj.*Culu£  Italia  uuLll  Luts  1st?  glim  lug  ui  Liie  2liu 
second,  when  the  process  continues.  The  scan  of  region  III  is  com¬ 
pleted  during  the  5th  second  and  starts  again  in  a  new  cycle  during 
the  6th  second.  The  same  occurs  for  region  II  (and  region  III  again) 
during  the  11th  and  12th  seconds,  respectively.  This  search  cycle  is 
shown  diagrammatically  in  Fig.  6.  The  order  of  search  within  a  fixed 
mode  is  to  increment  the  azimuth  steps  of  the  beam  position  (alter¬ 
nately  switching  between  antennas  #1  and  #2)  while  keeping  the 
elevation  fixed.  At  the  end  of  the  azimuth  scan,  the  elevation 
is  increased  by  one  step  and  the  azimuth  reset. 


COMPUTATION  OF  PROBABILITY  OF  DETECTION 


When  the  target  satisfies  the  criteria  for  being  inside  the  beam, 
i.e.,  6^  <  6t  £  bu  and  £  vT  £  v^,  a  probability  of  detection  Pdet 
is  computed  based  on  the  mode  of  the  radar,  the  target  range,  and  the 
position  of  the  target  relative  to  the  center  of  the  beam.  Preparatory 
to  this  computation  a  value  for  the  probability  of  a  false  alarm,  P^, 
must  be  selected.  This  is  the  probability  that  the  radar  system  will 
indicate  the  presence  of  a  target  at  a  particular  location  when,  in 

fact,  no  target  is  there.  For  this  model,  P^  has  been  set  equal  to 

-6  -9 

10  ,  although  it  will  probably  be  lowered  to  10  for  subsequent 

tests.  We  next  define  x  to  be  the  solution  of  the  equation 


P 


fa 


1 

m 


dt . 
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If  P,  "  10  ^  then  x  ■  4.7615  Is  the  solution.  We  also  define 

is 

R^  =  R^, /x,  where  j  is  the  mode  number  (1  ■  I,  II,  or  III) 

0  1 

We  now  suppose  the  target  is  inside  the  radar  beam  and  has 
coordinates  R^,  6^,  and  v^,  and  the  radar  is  searching  in  mode  j. 

We  compute  the  deviation  of  the  target  from  the  center  of  the 
beam  by  the  equation 

CT  ■  cos'1[sin  flT  sin  6c  +  cos  5T  cos  6c cos(vT-Vc) ] , 
where  6C=  ^(6^+fl^)  ®nd  vc  3  fe(v^-h>u)  are  the  elevation  and  azimuth, 
respectively,  of  the  beam  center.  The  effects  of  the  target  range 
and  deviation  from  beam  center  appear  in  the  value  of  the  parameter 
d,  which  is  calculated  as 


d  "  J^£  cot  (52.941€t). 

The  factor  52.941  insures  that  the  cosine  term  equals  — —  when 

o  ^ 

*  0.85  .  The  probability  of  detection  is  evaluated  as 


p,  -  JL  J 

det  /Tir  x-d 


.t*/2 
e  C  1  dt, 


The  value  of  Pdst  is  computed  by  using  an  approximation  for  the 
error  function  0  defined  by 


2  ry  -uz 

<t(y)  ®  -  e  u  du,  y  is  0. 

Jf  o 
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The  approximation  we  used  for  ®  is  found  In  [2,  p.  167]  and  la  the 
following: 

2 

$(y)  so  1  -  (a^  +  a2n2  +  a^3)  —  e“y  ,  y  z  0, 

/IT 

where  t]  ■  — - —  , 

i+py 

and  p  -  0.47047,  aL  -  0.3084284,  a2  -  -0.0849713, 


a3  -  0.6627698. 
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APPENDIX.  COMPUTER  TEST’S  OF  SIMULATOR 

Table  1  Indicates  the  results  obtained  from  20  computer  runs  of 
the  general  simulator,  consisting  of  four  sets  with  five  trial  runs 
in  each  set,  A  set  corresponds  to  one  of  four  possible  combinations 
of  target  height  h  and  standard  deviation  (see  definition  p.  8). 
These  latter  parameters  assume  the  values  h  ■  1.26  and  5,05  miles 
(target  crosses  horieon  at  ranges  of  100  and  200  miles,  respectively), 
and  o  »  0°  (target  heads  directly  for  projected  observer)  and  30° 
(target  heads  for  projected  observer  via  evasion  course).  For  all 
20  runs,  mean  free  time  t^  (see  definition  p,  7)  is  20  seconds  and 
target  speed  is  2000  mph. 

The  parameters  tabulated  for  each  run  are  (1)  the  time  that 
elapses  between  the  horiaon  crossing  of  the  target  and  its  detection, 
(2)  the  target  range  at  detection,  and  (3)  the  search  mode  in  which 
detection  occurred.  The  numerical  results  of  each  set  of  5  runs  are 
averaged . 
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Fig.  2  -  Representation  of  target-simulator  parameters  on 
surface  of  target  sphere 
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Fig.  3  -  Meridian  cross-section  of  earth  containing  target 
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